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ABSTRACT 


Title of Thesis: Predicting the Ignition Time and Burning Rate of 
Thermoplastics in the Cone Calorimeter. 


Name of degree candidate: | Donald Hopkins, Jr. 


Degree and Year: Master of Science in Fire Protection Engineering, 1995 
Thesis directed by: Dr. James G. Quintiere, Professor, Department of Fire 
Protection Engineering 


Ignition and burning rate data are developed for Nylon 6/6, Polyethylene, and 
Polypropylene in a Cone Calorimeter heating assembly. The objective is to examine a 
testing protocol that leads to the prediction of ignition and burning rate for thermoplastics 
from Cone data. The flame heat flux is not measured, but is inferred from Cone data. 
The constancy of the flame heat flux for thermoplastics in the Cone calorimeter is due to 
the geometry of the flame. The burning rate model is shown to yield good accuracy in 
comparison to measured transient values. 

Ignition and burning rate data are developed for Redwood and Red Oak in a Cone 
Calorimeter heating assembly. Measurements of the flame plus external heat flux are 
presented. The data is intended to be used for future work to develop a testing protocol 


and burning rate model for charring materials. 
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CHAPTER 1 INTRODUCTION 


Predicting the way a fire will behave under realistic conditions can be a 
challenging task. In order to examine the effects of fire on it’s surroundings and the 
hazard posed to occupants, information about the heat release rate, total energy, and 
flame spread characteristics of the fire need to be known. Subsequently, it is desirable to 
have a means of determining fire growth and spread in terms of measurable material 
properties. A general model for predicting the burning rate of materials is needed to 
accomplish this. 

Quintiere [1] developed a one-dimensional model which includes charring, 
vaporization, extinction, flame and heat conduction effects. However, unsteady 
solutions for the burning rate were not determined. Quintiere and Iqbal [2] developed a 
model that solves the one-dimensional unsteady heat transfer equations during the pre- 
heating and gasification periods using an integral method. They assume a polynomial 
profile for the temperature within the solid to satisfy the heat transfer boundary 
conditions. However, the effects of flaming are not included in the solution. 

The objective of this research is to develop transient burning rate models which 
utilize data obtained from the Cone Calorimeter [3]. The models will be dependent on 
the class of material; namely, thermoplastic, charring, dripping, and laminated. It is 
intended to first succeed for thermoplastic-like materials; materials approximated by 
constant surface temperature vaporization, which are large enough to be considered one 
dimensional in behavior. 

The ultimate goal of the research is to incorporate the burning rate models into 
fire growth simulations [e.g. 4, 5, 6, 7] to provide more pccuinie predictions and 
assessment of hazard. It is first necessary to succeed at developing a technique which 


utilizes Cone Calorimeter data to predict the transient burning rate, at least for thick 


thermoplastics and charring materials. The technique presented herein has been 
examined by Quintiere and Rhodes [8] for Black Polycast Polymethylmethacrylate 
(PMMA). The model is based on the formulation outlined by Quintiere [1] and 
implemented by Quintiere and Iqbal [2] for non-flaming pyrolysis of a thermoplastic. 
Further validation of the technique using Nylon 6/6, Polyethylene, and Polypropylene is 
attempted in the present study. The primary objective is to see if the method used by 
Quintiere and Rhodes [8] for PMMA is general for thermoplastic-like materials burning 
in the Cone Calorimeter. Furthermore, data is presented for charring materials; namely, 
Redwood and Red Oak to present insight into the extension-of the current model to deal 


with charring materials. 


CHAPTER 2 OVERVIEW OF THESIS 


2.1 Introduction 
The Cone Calorimeter is a commonly used device utilized to measure the mass | 


loss rate per unit area (mh") and the heat release rate per unit area (Q”) for a given 


constant external radiative heat flux. The ratio of the mass loss rate and the heat release 
rate gives the instantaneous heat of combustion (AH,) relative to the gaseous fuel 
produced during flaming combustion. It has been shown that the heat of combustion is 
generally constant for a material undergoing flaming combustion. The heat of 
combustion is measurable and is not expected to be scale dependent, i.e. change as a 
larger area of fuel is burned. Similarly the thermal and chemical properties of a 
decomposing material are also independent of scale, with the exception of multi- 
dimensional effects involving seams, joints, etc.. Conversely, the heat flux per unit area (q") 
to a materials surface is dependent on scale and on test conditions in the Cone 
Calorimeter. Subsequently, it can be shown that under one dimensional burning 
conditions, 


m" = f(q", thermo— chemical properties) (2.1) 


and, Op= miyAle (2.2) 


Equation (2.1) symbolically represents a model involving the surface heat flux and the 
needed properties. A specific model for a vaporizing non-charring thermoplastic model 
is examined. Experimentally, three materials are examined, namely; Nylon 6/6, 
Polyethylene, and Polypropylene. In addition, Redwood and Red Oak are examined 
experimentally to lay the foundation for a model to examine the behavior of charring 


materials. 


2.2 Burning Rate Experiments 

The burning rate experiments were conducted to measure the mass loss rate, 
ignition time, and surface temperature for both the thermoplastic and wood materials 
tested. A cone heater assembly with a continuous mass loss measurement was used for 
all of the experiments. Experiments were conducted until steady burning was reached. 
In addition, measurements of the flame plus external heat flux to the surface was 
measured for the wood samples using a 3 mm diameter heat flux gage. This was not 
done for the thermoplastic materials. However, previously Rhodes [9] measured the 
surface heat flux for PMMA in the Cone Calorimeter and concluded that the flame heat 
flux (convection plus radiation) 1s constant for a thermoplastic-like material burning in 
the Cone Calorimeter. The flame heat flux can be deduced from the steady burning rate 


data for the thermoplastics in which Equation (2.1) becomes, 


ao 
mh" => (2.3) 
where, q” is the net surface heat flux, 
and, L is the effective heat of gasification. 


At steady burning, the incident flame heat flux is deduced from the net surface heat flux 
in equation (2.3). The gage measurement for the thermoplastic materials would have 
only provided a second means of determining the flame plus external heat flux to the 
material surface. The current results verify that the flame heat flux for a thermoplastic- 
like material burning in the Cone Calorimeter is approximately constant. However, due 
to charring and the associated change in flame height for wood burning in the Cone 
Calorimeter, the constancy of the flame heat flux does not appear to be one of the 


attnibutes of the burning behavior. 


2.3 Application of Models 

The data obtained from the experiments was examined in conjunction with 
ignition and burning rate models designed to yield the needed thermo-chemical 
properties, i.e. the thermal conductivity, density, specific heat, and heat of gasification. 
The properties are model dependent, but are not void of physical significance. 

The experimental and modeling details will be reviewed and their results 
presented. A recommended testing protocol for thermoplastics is presented in Chapter 5. 
The method used to elucidate the needed properties for the wood materials is applied, but 
currently a model has not been used to attempt to predict the ignition time and burning 
rate for wood. The results of the wood experiments are presented for aid in the future 
development of a model for charring materials, which will be based on a solid phase 
model developed by Quintiere [1]. 


CHAPTER 3 THERMOPLASTIC EXPERIMENTAL SET-UP, 
PROCEDURE, AND MEASUREMENTS 


3.1 Experimental System 

| Ignition and burning rate experiments for the thermoplastic materials were 
performed using a radiant cone heater assembly. The apparatus, shown in Figure 3.1 
[9], consisted of a cone heater, a load cell, a methane pilot ignitor, and a data acquisition 
system. A scanner and a voltmeter were used to monitor thermocouple, heat flux 
transducer, and load cell measurements. A Hewlet Packard (HP) computer was utilized 
as the data acquisition system. As data was measured it was transferred to a second 
computer where it was recorded. The voltage measurements corresponding to 
temperature, mass, and heat flux were converted to the appropriate metric units prior to 
their recording. The HP Basic program used for the experimentation is included in 
Appendix A, although a description of the program is not included in this report. 


3.2 Thermoplastic Samples 

Three thermoplastic-like materials were selected for evaluation in this study; 
namely, Nylon 6/6, Polyethylene, and Polypropylene. These materials were selected 
based on availability from commercial retailers. The samples were nominally 10 cm (4 
in.) x 10 cm (4 in.) x 2.5 cm (1 in.) thick. Previously, experiments were carried out by 
Rhodes [9] using Black Cast Polymethylmethacrylate (PMMA). The current testing and 
analysis is an attempt to generalize the previous results found for PMMA. 

Samples were placed on a standard cone metal holder in the horizontal orientation 
on a bed of Kaowool. The Kaowool was used to insulate the back side of the specimen 
to minimize heat loss effects. In order to maintain only one dimensional burning, 


cardboard was bonded to the sides of the samples to inhibit edge burning effects. 
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Figure 3.1 Schematic layout of burning rate apparatus. 


Furthermore, aluminum foil was wrapped around the edges and the back of the sample to 


prevent dripping. 


3.3 Experimental Procedure 

The experimental procedure consisted of exposing a sample, in the horizontal 
orientation, to a constant external irradiance from the cone heater assembly. The time to 
piloted ignition was measured and mass loss data recorded for each test. Prior to the 
onset of each experiment the cone heater needed to be set to the desired constant external 
irradiance. The initial incident heat flux was determined using a Medtherm 1 inch 
thermopile type heat flux transducer (Model # 64-5SB-20) situated such that it was in 
the same location that the center of the sample would be, and so that the face of the gage 
was 1 inch below the base of the cone heater. The experiment was not started until a 
constant heat flux recording was obtained for at least a one minute period. The heat flux 
gage could not be left in place during the experiment. Subsequently, the constancy of the 
cone heater was ensured by monitoring the temperature of the cone coil. Two 
thermocouples were located in the cone heater coil. The thermocouple readings were 
averaged and displayed every two seconds so that the cone heat flux, which corresponds 
a given coil temperature, could be monitored. The cone coil temperature was kept 
constant (+ 5° C) by manually adjusting the current to the coil using a 220 Volt 
transformer. This approach allowed for adjustment of the current to compensate for 
temperature fluctuations induced by flame impingement on the cone coil. Figure 3.2 


shows a typical cone temperature variation for Nylon subjected to a 75 kW/mz2 exposure. 
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Figure 3.2 Cone heater coil temperature variation with a 75 kW/m? external irradiance. 


3.3.1 Ignition 


The arrangement of the assembly was such that the top of the sample was initially 
located 1 inch below the base of the cone heater as is done for standard Cone Calorimeter 
tests [3]. The load cell was oriented between two guide bars to ensure proper placement 
of the sample underneath the center of the cone heater assembly. A one inch methane 
flame, located on one edge approximately 1/2 inch above the surface of the sample, was 
used as a pilot ignition source. Figure 3.3 shows the arrangement of the pilot and the 
assembly. 

Once the cone was set to the desired constant external irradiance and the sample 
was prepared and situated on the load cell, the load cell was shifted into position 
underneath the cone. Simultaneously, the data acquisition system was initiated. It should 
be noted that sliding the sample into place caused some bouncing of the load cell. This 
effects the mass loss recordings for the first few seconds of the experiment. This is not 
however expected to have adverse effects on the results since ignition nominally took 
longer that ten seconds. 

The ignition time is defined as the time at which a continuous flame is supported 
on the material surface. In some instances flashing occurred on the surface of the sample 
prior to sustained flaming. However, in all cases the ignition time was taken as the time 
at which flaming was sustained over the entire surface of the specimen. Ignition times 
for all of the thermoplastic tests for irradiances ranging from 20 kW/m2 to about 80 
kW/m2 are shown in Appendix B. 

In some experiments, a fine wire Type K thermocouple (0.003 in. diameter) was 
mounted on the thermoplastic material surface to attempt to measure the surface 
temperature at ignition. The thermocouple wire was placed on the sample surface by 
heating the wire so that it recessed into the surface of the material. Measurements of the 


surface temperature were recorded until the onset of burning. Figure 3.4 shows a typical 
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Figure 3.3 Experimental burning rate apparatus. 
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Figure 3.4 Surface temperature of Polyethylene as a function of time with a 36 kW/m2 
external irradiance. 
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result of the measured surface temperature for the preheating period for Polyethylene 
exposed to a 36 kW/m? external irradiance. The temperature at the onset of the sudden 
rise in temperature is defined as the ignition temperature (Tj,). Results for the surface 
temperature measurements for the thermoplastic experiments at various external 


irradiances are shown in Appendix C. 


3.3.2 Mass Loss Rate 

Experiments to determine the mass loss rate of the materials were performed 
concurrently with the ignition experiments. The mass loss of the samples was recorded 
using a load cell. The mass loss readings were recorded every two seconds for 
approximately a 1200 second period. 

The transient mass loss rate is found using a five point least squares fit of the 


mass loss data. The mass loss rate at some time t; is given as, 


i+2 i+ 2 i+ 2 
ect), Cm) ae) 
f="I—2 n= i-2 n= i-2 
m" — 
i+ 2 a) i+2 2 
Bae erste) | y ) 
n= i-2 N=i-2 
where, i is a given measured data point, 
m is the mass at 1, 
and, t is the time at 1. 


Figure 3.5 shows a typical mass loss rate result for Nylon exposed to an 80 kW/m2 
irradiance. Results for Nylon at other exposures, Polyethylene, and Polypropylene are 


shown in Appendixes D, E, and F, respectively. 
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Figure 3.5 Transient mass loss rate of Nylon with an 80 kW/m2 external irradiance. 
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3.4 Calibration 


Calibration of the load cell, thermocouples, and 2.5 cm (1 inch) heat flux gage 
was necessary before any experiments were conducted and only occasionally once 
testing began. The methods used to calibrate each of the instruments are discussed 


below. 


4.1 Load Cell 
The load cell was calibrated by first adjusting the weight of the sample holder and 

insulation to approximately zero. At this point, standard weights between 1 and 100 
gram were added to the load cell and the voltage readings were recorded. The 
relationship between weight and voltage was found to be linear, with a calibration 
constant of 582.0 g/V. Because only mass loss rate is desired, only the difference in 
mass is necessary. Therefore, small fluctuations in the zero weight did not affect the 
results. At the start of each day, the load cell was checked for accuracy using 1 and 2 
gram weights. Furthermore, in some instances the mass of the sample was measured 
using a scale and the results compared to that indicated by the load cell. No significant 


errors were found for the mass loss measurements. 


3.4.2 Heat Flux Gauge 
The 2.5 cm (1 inch) heat flux gauge was initially calibrated using the standard 


NIST heat lamp calibration box. It was tested with water at room temperature and again 
with water at 46 °C. No notable differences in the calibration constant were observed. 
Although the water temperature did not affect the calibration constants, it did affect the 
zero reading. Unlike the mass loss, the heat flux measurement must be zero when there 
is no incident heat flux. The hot water resulted in an initial negative heat flux 


measurement due to convection between the sensor and the atmosphere. The heat flux 
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measurement was calibrated to be zero, under no external irradiance, before each 
experiment in order to obtain consistent results. Besides the zero flux measurement, the 
1 inch gage did not require calibration before each of the experiments because it was 
positioned such that it did not get disturbed or touched. The calibration constant for the 1 


inch heat flux meter was determined to be 5235 kW/m2 V. 


.4.3_ Thermocouples 

The thermocouples used to measure the coil temperature and surface temperature 
were calibrated before any testing began. This was done by placing them in an ice bath 
and also in boiling water for a prolonged period of ttme. Temperature readings were 
found to be accurate within +2 °C. The thermocouples in the coil did not require 
additional calibration as the absolute temperature was not necessary. The thermocouples 
for the surface were calibrated against each other before each experiment. 

Although the experiments yielded extraneous results for some tests, no major 
calibration errors were found. Inconsistencies in the results are likely to be due to factors 


other than calibration errors. 
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CHAPTER 4 MODEL FOR THERMOPLASTICS 


The model used for this analysis is based on the formulation of Quintiere [1] as 
presented by Quintiere and Iqbal [2]. The model is one-dimensional and assumes that 
surface vaporization occurs at a specified temperature, Ty. The model can be extended to 
include an analysis for charring materials, but at the present time only thermoplastic 


(non-charring) materials will be considered. The model will be briefly described below. 


4.1 Preheating to Ignition 
The one-dimensional unsteady heat conduction equation applies to the preheating 


period. Constant properties are assumed such that the governing equation is, 


oT Cal 

anes a ay? (4.1) 
with a constant initial temperature , 

T=Tp@t=0. (4.2) 


Furthermore, considering convection at the surface (y = 0) and radiative heat loss, 


oT (4.3) 


—kZy | eo 4" = €4’e — B, T-To) - eoT” 


An approximate integral solution is applied to formulate a solution to the above problem. 


Therefore, the solution to Equation (4.1) comes from integrating the equation between 


zero and some penetration depth, 5 , where at y = 6, 


and, — = 0. (4.4) 
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A quadratic profile is assumed for T such that the three conditions given by equations 
(4.3) and (4.4) are satisfied: 


, 2 
itscyi be % (!-3) Badin 3 


Applying the above boundary condition and integrating equation (4.1) from 0 to 6 yields, 
d 2 
< (q°8") = 60" (4.6) 


If q" is assumed to be constant, which is a reasonable assumption for large external 


irradiances, then equation (4.6) yields, 
5 = aj 6cxt (4.7) 


A more complete solution to the above problem has been shown by Abu-Zaid and Atreya 


[10] to yield; 


5= (4.8) 


where e€4 was shown to vary between 1.6 at 15 kW/m2 and 1.9 at50 kW/m2. Taking 


the asymptote as 2.0 yields equation (4.7). Since the solution presented here is 
approximate, the accuracy sacrificed by using equation (4.7) is assumed to be acceptable, 
although the error will be greatest, about 10%, at low external mrradiances. 


Substituting the equation for 5 from equation (4.7) into equation (4.5) at y =0 


produces, 


"5 A 
peengors este th (4.9) 
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Solving equation (4.9) for tis gives, 


(4.10) 


| This approach also gives a method for determining the surface temperature as a 
function of time. This is done implicitly by selecting T,, calculating the corresponding 
net heat flux, and using equation (4.10) to determine the time. This result also allows for 
the determination of a critical flux for ignition, q" «, by extrapolating ignition data for 


(tig)-1/2 to zero. At this intercept, 


1 
oh aa th [h.(T,,- T») a eoT* | = q'. (4.1 1) 
E 


The critical flux for ignition is determined using equations (4.4) and (4.10). Once the 
critical flux is found, the temperature at the critical flux can be determined using equation 


(4.11). 


4.2 Burning Rate 

The governing equations for the gasification period can be derived in integral 
form by selecting the appropriate control volumes (CV; and CV>2) for the vaporization 
plane and the solid, as shown in Figure 4.1 [2]. The governing equations for the 
burning rate follow from equation (4.1) which governs conduction to the material below 


the vaporizing plane at a fixed temperature. Subsequently, at y = 0, 
T=T, (4.12) 
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fixed reference 


;, gas, y=0 


Ck OT/dy) 4 y=0 


d(y+ 3)/dt 


solid is at rest 


Figure 4.1 Heat and mass transfer processes for a thermoplastic-like solid fuel. CV1 is 
around the vaporizing interface. CV2 is bounded by the interface and 


thermal penetration depth, 6. 
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and, 


J) T " o. tf 
= See = q' — m"AH, 
where, om" is the mass loss rate per unit area, 


AH, is the heat of vaporization, 


w 


and q is the net surface heat flux. 


It should be noted that the net surface heat flux, q", is different for the gasification 


(4.13) 


period in equation (4.13) than for the preheating period as defined in equation (4.4). 


The net surface heat flux for the gasification period is, 
q"= eat 4" y~ eT, 


v 


where q" g is the flame heat flux, which is defined as, 


q".= eq" ..+ G' ae 


(4.14) 


(4.15) 


Where, q",, is the flame incident radiant heat flux, and q",..is the convective flux, 


from a stagnation film model [1], and can be shown to be, 


7 aoe cul T 4.16 
Vass ate) | Yokel Te) (4.16) 
m'c, 
where, a > 
( & ) is the mass transfer ’blocking factor” that is 1 for m" ->0, 
Cr= 
v6 is the ambient mass fraction of oxygen, 
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re is the flame radiative fraction, 


AH, is the heat of combustion, 
iy . is the stoichiometric fuel to oxygen mass ratio, 
Cg is the gas phase specific heat at constant pressure. 


An alternate way of expressing this is, 


q"ac= h.(—2—)(T,-T,) (4.17) 
e&-1 


where Tg is an effective flame temperature. 
Assuming a quadratic profile for the temperature, which satisfies the boundary 
conditions supplied by equations (4.2) and (4.12), 


T—2 
Ten 


v 


y AD 
= (l-— 4.18 
( 5 (4.18) 


and substituting the profile into equation (4.13) yields, 


m" AH, =q" = (leestths (4.19) 


where the last term represents transient conduction heat loss into the solid. Integrating 


equation (4.1) from y = 0 to y=4, as was done for the ignition analysis, and 


substituting in the profile from equation (4.18) yields, 


Sens ck Sere - (4.20) 


where at ignition, T = T, and 6 = | 6at = bis. (4.20a) 
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If the total flame heat flux is assumed to be constant, or if the dependence of the buming 


rate on the blocking factor is ignored, then equation (4.20) can be solved exactly. It can 


be shown that, 
2 
6. AH, 6,.—5 w 6-95 
t-t, = —-——| 2 - (4.21) 
* 6a - | 8 5,-5,, 
where, 6. = aut a Steady value, (4.22) 
and, L = AH, + c(T, — T>), is the heat of gasification. (4.23) 
It follows that the steady state mass loss rate is given by, 
rh"= (4.24) 


It should be noted that in addition to the flame heat flux, other properties are 
needed to obtain a solution. These properties need to be derivable in a convenient 


manner consistent with the burning rate and ignition models. The properties are, 


€,p,c, k, T, or T;,, and AH, or L. 


4.3 Flame Heat Flux 

For the present experiments, only the steady mass loss rate measurements are 
used to determine the flame heat flux. Rhodes [9] used surface heat flux measurements 
in addition to the steady mass loss measurements for his study with PMMA. In 
addition, a simulated sample was used to examine the effect of the flame heat flux. This 


will be discussed. 
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4.3.1 Flame Radiation 


The emissivity of a flame can be represented as, 


e=1- cane (4.25) 


where, K is the absorption coefficient, 


and _1,, is the mean beam length. 


According to Orloff and deRis [11] , k and 1, can be computed from an algorithm for 


pool fires. For cases where the flame height 1s greater than twice the diameter of the 
base, which is the case for all of the thermoplastic-like materials tested, the mean beam 
length is, 

]_ = 0.65D (4.26) 
where D is the diameter of the base, which is 10 cm for the materials tested in the cone 
heater assembly. This says that the mean beam length is constant for a given diameter 
when the flame height is greater than twice the diameter. It follows that the flame 
emissivity is constant and subsequently the radiative flame heat flux is approximately | 


constant for a constant flame temperature because, 
q"a= €,0T. (4.27) 
In addition, the transmissivity is defined as, 
Tze," (4.28) 
Because the mean beam length is constant, it follows that the transmissivity is constant. 
This suggests that the amount of radiation from the cone heater that reaches the surface is 


constant for a thermoplastic material. Rhodes [9] found that the transmissivity of PMMA 


is approximately 0.91, suggesting that the flame is very transparent. 
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4.3.2 Convective Flame Heat Flux 
The convective heat flux, assuming a blocking factor of 1, can be estimated from 


equation (4.16) considering, 


ho = 10 W/m2K AH /r = 13 kI/g 
cp =1.0J/gK x, = 0.4 
Y,.. = 0.233 


oxX@ 
A maximum q" f.¢ = 15 kW/m? is determined since blocking would reduce this value as 


the burning rate increases. 


4.3.3 Total Flame Heat Flux 
Assuming the surface emissivity , €, 1s approximately 1, equation (4.15) 


indicates that the total flame heat flux is approximately constant. This constancy has 
been shown for PMMA by Rhodes [9] and has good utility in aiding in the analysis of 
Cone Calorimeter data. 

Measurements of the flame plus external heat flux were not attempted for the 
thermoplastics examined due to difficulties encountered in the PMMA experiments. 
However, it was shown by Quintiere and eines [8] that the flame heat flux for 
thermoplastic-like materials burning in the Cone Calorimeter is constant. This was 
accomplished through measurements of the flame plus external heat flux to the material 
surface during the burning rate experiments and by utilizing a simulated sample methane 
gas burner. Since the phenomenon of the constant flame heat flux is critical to the 
analysis presented within for thermoplastics, a brief overview of the procedure using the 
simulated sample is presented here. 

A simulated sample was used to examine the effect of flame heat flux on burning 


rate. The simulated sample was steel and measured 10 cm x 10 cm x 7.5 cm thick. A 


2. 


14.3 mm (9/16 inch) hole was located in the center of the sample to allow for heat flux 
measurements using a 3.175 mm (1/8 inch) heat flux meter. The burner was filled with 
glass beads and erate was injected through the base. The glass beads insured a 
steady flow rate of methane through the top surface, with the rate of flow of the gas 
being controlled by a calibrated flow meter and a regulator. 

The methane burner simulated sample was positioned beneath the cone heater in 
exactly the same position as the thermoplastic samples. The sample was subjected to a 
50 kW/m2 constant external irradiance. Considering that the sample was utilized to 
examine the effects of flame heat flux in the Cone Calorimeter, the sample was placed 
under the cone, not in the open environment, since the assembly tends to induce a 
specific airflow which influences the burning rate. The sample was exposed to a 
constant external flux so that flame impingement on the cone coils did not change the heat 
flux to the surface. As was previously discussed, the cone coil temperature was 
monitored and could be adjusted using a transformer. Figure 4.2 shows the steady state 
sensor heat flux measurements for an external exposure of 50 kW/m2. The results show 
a constant flame heat flux above the external radiant heat flux. This is approximately 27 
kW/m? for the methane over a range of energy release rates per unit area of 200 - 600 
kW/m2, which was determined from the methane flow rate and the heat of combustion of 
methane (50 kJ/g). 

The constancy of the flame heat flux may be attributed to the shape of the flame 
for materials burning in the Cone Calorimeter. Figure 4.3 shows the flame for PMMA 
burning in the Cone. The long column shape of the flame is typical of all of the 
thermoplastic-like materials tested. The flame height to effective diameter is on the order 
of 4. 
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Net Surface Flux of Methane Burner 
exposed to a 50 kW/m‘2 External Flux 


Flux (kW/m42) 
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Heat Release Rate (kW/m42) 


‘Figure 4.2 Steady state heat flux, measured by sensor, of methane gas burner with a 50 
kW/m? external irradiance. 
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Figure 4.3 
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CHAPTER 5 ANALYSIS OF THERMOPLASTIC RESULTS 


The analysis of the results is done in a manner consistent with the protocol 
outlined by Quintiere and Rhodes [8] for PMMA. The intent is to see if the protocol is 
general for all thermoplastic-like materials. Results obtained from the thermoplastic 
experiments are used to obtain the required properties to utilize the ignition and burning 
rate models. While most of the necessary properties are deduced from the experimental 
data, the density is determined by measurement, and the thermal diffusivity, which is 
assumed to be approximately constant for a given material, is obtained from the literature 
[12, 13]. The density and the thermal diffusivity for the three thermoplastics are shown 
in Table 5.1. 

The approach used to obtain the needed properties will be demonstrated for 
Nylon 6/6. The approach is identical for the analysis of the other materials. The 


properties for all of the materials can be seen in Table 5.3. 


Table 5.1 


Measured and Literature Properties 


Polyethylene aS elec | 
| eee 10-7 aa 10-7 | 6.736 x 10-8 
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5.1 Ignition 

In order to predict the ignition time and the critical flux for ignition, the ignition 
temperature and the thermal inertia must be determined. Once this is accomplished it is 
possible to determine the critical flux and the ignition time as a function of the external 


irradiance. 


5.1.1 Ignition Temperature and Thermal Inertia 


Figures 5.1 - 5.4 show plots of external irradiance as a function of ignition time 
for Nylon 6/6, Polyethylene, Polypropylene, and PMMA, respectively. The line 
through the data represents the best fit of the data below 40 kW/m2 where the results 
should be better since for lower heat fluxes ignition takes longer. The best fit line in 


Figure 5.1 gives an intercept of approximately 14 kW/m2, which represents the critical 


flux for ignition for Nylon 6/6. The values for the critical flux for ignition, q_ , for the 


other materials are shown in Table 5.3. Using the value for the critical flux, equation 


(4.11) can be used to determine the ignition temperature, Tj,, given To = 20 °C, h, = 10 
W/m2-K (determined from an analysis for natural convection over a flat plate), and € = 


1.0. This yields an ignition temperature for Nylon 6/6 of T;, = 380 °C. While this value 
is lower than the measured ignition temperature for Nylon, it will be shown that the 
deduced value yields good results when used in the model. Table 5.2 shows a 
comparison of the measured and calculated ignition temperatures for the materials tested. 
For all of the materials, the measured ignition temperatures were higher than the deduced 
values. While this may be an attribute of the modeling procedure, there is also 
uncertainty in the measured values, especially for Nylon where only one result was 
obtained. 


Using the slope of the line in Figure 5.1, along with equations (4.3) and (4.10) 
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External Flux (kW/m‘2) 


Figure 5.1 Ignition data for Nylon. 
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Figure 5.2 Ignition data for Polyethylene. 
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Figure 5.3 Ignition data for Polypropylene. 
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Figure 5.4 Ignition data for PMMA. 
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80 


the thermal inertia can be calculated as follows: 

From equations (4.3) and (4.10 ), 

h.(T,,—T,) +£0T;, 
2 
$kpe (T;,—T») 


q ex C551) 


Sri € 
ig 2 
3 kpe (T,,—To) 


Subsequently, from equation (5.1), the slope of the line in Figure (5.1) is, 


E 
Slope = on (S22) 
3 kpc (T;,— Ty) 
Equation (5.2) can be solved to determine the thermal inertia, 
iz 
sepdres walt aAtesee ete nst 5.3 
Pe = 2 | Slope(T,,-T) Se 


Therefore, for Nylon 6/6 , the slope from Figure (5.1) is 0.00364 and subsequently the 
thermal inertia is, 
kJ’ 
kpc = 0.874 ——; 
fe m‘s K? 


The values for the thermal inertia, kpc, for the other materials are shown in Table 5.3. 


Figures 5.5 - 5.8 show plots of ignition time versus external flux for Nylon, 
Polyethylene, Polypropylene, and PMMA, respectively, along with the results of the 
model using the deduced values for the ignition temperature and the thermal inertia as 


input for equation (4.10). 
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5.1.2 Thermal Conductivity and Specific Heat 

Knowing the thermal inertia, kpc, as was determined above, the density, p, and 
the thermal diffusivity, « = k/pc , which varies only slightly with temperature, the 
thermal conductivity, k, and the specific heat, c, can be determined. 
For Nylon 6/6 with p = 1169 kg/m3 and a = 1.24 x 10-7 m2/s, the thermal conductivity 


and the specific heat are, 
k = 3.29 x 104 kW/m K 
and, c= 2.27 kI/kg K 


Although the thermal diffusivity must be independently determined, the above method 
gives a means of determining the specific heat and the thermal conductivity. Table 5.3 
shows the specific heat and the thermal conductivity for the materials tested. 


Table 5.2 
Calculated and Measured Ignition Temperatures . 


[ Material ‘Calculated Ignition Ignition | Measured I Ignition 
— (°C) Temperature ture((C) C) 
PC Rne 6/6 96 eae a = 500 
250 - 355 
Dpeimeteomemnihcemieieeiye 
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Figure 5.5 Calculated and measured ignition time of Nylon as a function of external 
irradiance. 
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Figure 5.6 Calculated and measured ignition time of Polyethylene as a function of 
external irradiance. 
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Polypropylene 
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Figure 5.7 Calculated and measured ignition time of Polypropylene as a function of 
external irradiance. 
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Figure 5.8 Calculated and measured ignition time of PMMA as a function of external 
irradiance. 
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5.2 Burning Rate 

Since most thermoplastic-like materials approximate vaporizing solids, it is 
possible to represent the steady-state mass loss data using equation (4.24). Therefore, if 
the flame heat flux is assumed to be constant, which has been shown to be the case for 
thermoplastic-like materials burning in the Cone Calorimeter [8], then equation (4.24) 


can be written as, 


1 (q",—€0T,) 
m" = ()d" ex + ed (5.4) 


Subsequently, a plot of the steady state mass loss rate data as a function of the external 
flux has utility in determining the heat of gasification and the total flame heat flux. 
Figure 5.9 - 5.12 show plots of steady state mass loss rate versus external irradiance for 


Nylon 6/6, Polyethylene, and Polypropylene, and PMMA, respectively. 


5.2.1 Effective Heat of Gasification and Total Flame Heat Flux 
Equation (5.4) suggests that the slope of the fit to the data in Figures 5.7 - 5.9 


represents the inverse of the heat of gasification, 1/L. Subsequently, for Nylon 6/6, 


ins 5. 
(+-) = 0.26456 +5 


Therefore, 


L = 3.78 kJ/g. 
Furthermore, the intercept of the fit to the data in Figure 5.7 is, 


q" ,—e0T, 


v 


ie me5.5) 


For Nylon 6/6 this is, 
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Therefore, for an average vaporization temperature, taken to be the deduced value of the 
ignition temperature (380 °C), and an emissivity, € , assumed to be 1.0, the flame heat 
flux is estimated using equation (5.5). For Nylon the flame heat flux is, 


q's $7), 28 
m? 


Table 5.3 shows the heat of gasification and the flame heat flux for the other materials. 


Table 5.3 
Deduced Properties 


i 3.29x 104 | 639x104 | 3.81x104 |0.432x103| kW/mK 


kJ/kg K 


Pap 3780 3580 3070 2770 


qa TF kW/m? 
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Figure 5.9 Steady state mass loss rate of Nylon as a function of external irradiance. 


43 


Polyethylene 


Slope = 0.27916 g/kJ 


Mass Loss Rate (g/m“2-s) 


0 10°20 9 SU reee Ue" OU tO Usb .\ OU ae oU 


Heat Flux (kW/m42) 


Figure 5.10 Steady state mass loss rate of Polyethylene as a function of external 
irradiance. 
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Figure 5.11 Steady state mass loss rate of Polypropylene as a function of external 
irradiance. 
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Figure 5.12 Steady state mass loss rate of PMMA as a function of external irradiance. 
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5.2.2 Transient Mass Loss Rate 


The mass loss rate per unit area can be calculated using the properties in Tables 
5.1 and 5.3 along with equations (4.19) and (4.21). The convective heat transfer 


coefficient, ho, was taken to be 0.01 kW/m2-K (determined from an analysis for natural 


convection over a flat plate). 

Figures 5.13 - 5.18 show the calculated and experimental mass loss rates for 
Nylon exposed to range of external irradiances from 31 - 80 kW/m2.. It should be noted 
that at low external heat fluxes there are discrepancies in the calculated versus the 
experimental ignition times for Nylon 6/6. This can be seen for the exposure to a 31 


kW/m2 irradiance, shown in Figure 5.13. The reason for the discrepancy is that at low 


heat fluxes (q".., < 35 kW/mz2) a bubble formed on the surface of the Nylon sample, 


thus preventing the release of pyrolisates. Ignition did not occur in these cases until the 
bubble over the surface broke. This phenomenon was only seen in the Nylon 
experiments and only at low external heat flux exposures. Once the Nylon samples 
ignited they quickly reached a steady state burning rate. While the Nylon was burning 
small bubbles from the surface could be seen bursting and burning in the air near the 
sample. With the exception of the erroneous ignition time shown in Figure 5.13, the 
calculated mass loss rate results for Nylon are in good agreement with the experimental 
data. 

The quick rise to a steady burning rate after ignition, as well as an approximately 
constant flame height, which is much taller than the top of the cone heater, is indicative 
of the burning behavior of all three of the thermoplastic materials tested. Polyethylene 
and Po vmronviene both seemed to become a pool of liquid inside of the aluminum foil 
while the samples were undergoing steady burning. In some instances the material could 
be seen boiling inside of the flame. 


Calculated and experimental mass loss rates for Polyethylene exposed to external 
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irradiances ranging from 36 - 87 kW/m? are shown in Figures 5.19 - 5.23. Fluctuations 
can be seen in some of the Polyethylene experimental data. The exact reason for the 
fluctuations is unknown, but it is assumed that they were caused by currents in the lab 
causing bouncing of the load cell. Despite the fluctuations, the calculated mass loss rate 
results for Polyethylene appear to be in good agreement with the experimental data. 

Calbitatea and experimental mass loss rates for Polypropylene exposed to 
external irradiances ranging from 27 - 65 kW/m? are shown in Figures 5.24 - 5.30. The 
calculated mass loss rate results appear to be in good agreement with the experimental 
data. 

Calculated and experimental mass loss rates for PMMA exposed to external 
irradiances of 25, 50, and 75 kW/m2 are shown in Figures 5.31 - 5.33, respectively. 

It should be noted that in Figures 5.17, 5.21, 5.22, 5.23, 5.28, and 5.30, the 
mass loss rate rises above the steady value after some duration. This rise in mass loss 
rate is caused by an increase in temperature of the material associated with the fact that 
there is only a fraction of the steady state penetration depth remaining and subsequently 
less material is being heated. Because the sample is insulated with Kaowool on the | 
backside, little or no heat is lost through the back, and consequently the mass loss rate 
increases. This phenomenon is expected to have been present in all of the experiments, 


prior to burnout, should the test have been allowed to proceed long enough. 
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Figure 5.13 Calculated and experimental transient mass loss rate of Nylon with a 31 
kW/m2 external irradiance. 
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Figure 5.14 Calculated and experimental transient mass loss rate of Nylon with a 35 
kW/m? external irradiance. 
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Figure 5.15 Calculated and experimental transient mass loss rate of Nylon with a 50 
kW/m2 external irradiance. 
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Figure 5.16 Calculated and experimental transient mass loss rate of Nylon with a 60 
kW/m2 external irradiance. 
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Figure 5.17 Calculated and experimental transient mass loss rate of Nylon with a 75 
kW/m? external irradiance. 
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Figure 5.18 Calculated and experimental transient mass loss rate of Nylon with an 80 
kW/m? external irradiance. 
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Figure 5.19 Calculated and experimental transient mass loss rate of Polyethylene with 
a 36 kW/m? external irradiance. 
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Figure 5.20 Calculated and experimental transient mass loss rate of Polyethylene with 


a 59 kW/m2 external irradiance. 
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| Figure 5.21 Calculated and experimental transient mass loss rate of Polyethylene with 
. a 61 kW/m? external irradiance. 7 
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Figure 5.22 Calculated and experimental transient mass loss rate of Polyethylene with 
a 70 kW/m? external irradiance. 
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Figure 5.23 Calculated and experimental transient mass loss rate of Polyethylene with 
an 87 kW/m? external irradiance. 
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Figure 5.24 Calculated and experimental transient mass loss rate of Polypropylene with 
a 20 kW/m? external irradiance. 
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Figure 5.25 Calculated and experimental transient mass loss rate of Polypropylene with 
a 27 kW/m? external irradiance. 
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Figure 5.26 Calculated and experimental transient mass loss rate of Polypropylene with 
a 34 kW/m2 external irradiance. 
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Figure 5.27 Calculated and experimental transient mass loss rate of Polypropylene with 
a 39 kW/m? external irradiance. 
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Figure 5.28 Calculated and experimental transient mass loss rate of Polypropylene with 
a 50 kW/m? external irradiance. 
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Figure 5.29 Calculated and experimental transient mass loss rate of Polypropylene with 
a 61 kW/m? external irradiance. 
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Figure 5.30 Calculated and experimental transient mass loss rate of Polypropylene with 
a 65 kW/m? external irradiance. 
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Figure 5.31 Calculated and experimental transient mass loss rate of PMMA with 
a 25 kW/m? external irradiance. 
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Figure 5.32 Calculated and experimental transient mass loss rate of PMMA with 
a 50 kW/m? external irradiance. 


68 


PMMA - 75 kW/m‘2 


a = A ot aed > PS - bee tee a= PP 
sels eae vicwe-e tt eM tsS oiany 
7. o 


eDC \jGec-0 
og: 2-pjet 


is 
< 

£ 
= 
2 
is°} 
~% 
o 
a 
Z 
= 

© Experiment 
= Calculated 
0 100 200 300 400 500 600 
Time (sec) 


Figure 5.33 Calculated and experimental transient mass loss rate of PMMA with 
a 75 kW/m? external irradiance. 
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CHAPTER 6 DIMENSIONLESS FORM OF SOLUTION 
FOR THERMOPLASTICS 


An alternate way of examining the thermoplastic solution is to look at it in terms 


of dimensionless parameters. The governing equation for the calculation of the ignition 


time and the burning rate can be represented in dimensionless form. The following 
dimensionless variables are defined: 


mh"L 
peut 6.1 
a | (6.1) 
AH, 
M) — L (6.2) 
8 
x=— (6.3) 
5, 
Using Equations (4.20a), (4.22), and (6.2) it can be shown that 
6, a, 
=# = fuse, = 0-9 fa (6.4) 


Using these dimensionless parameters, from equations (6.1) - (6.3), the mass loss rate 


per unit area, given by equation (4.19), is represented as, 


aie P = iN a 1 (6.5) 


The dimensionless time is given as, 


(6.10) 


6 is a parameter that physically represents the thermal penetration into a semi-infinite 


solid. 


The steady state value, from equation (4.22), is given as, 


Ole a (6.7) 


56. = 4 6at... (6.8) 


The time for ignition, from equation (4.10) is approximately given by 


GyrTe 
q” 


a Te ee (6.9) 


Subsequently, the dimensionless time after ignition , from equation (4.21) is, 


1-1, = 3 (4 1.5t,,.-x} — In fie (6.6) 


The net surface heat flux, q", changes at ignition to account for the flame heat 
flux, 


q",» which includes both radiative and convective fluxes to the surface. For our 


application, it has been shown that this net flame flux is constant for each thermoplastic 
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regardless of time and external irradiance [8]. The net surface heat flux before ignition is 


represented as, 


d= G" eq — hy (Tig- To) — OT;, for t< ti, (6.11) 


and after ignition as, 


4, = 4" ex + 4", — OT;, for t>t, (6.12) 


Figure 6.1 shows a plot of dimensionless mass loss rate versus dimensionless 


time after ignition for a constant value of external heat flux over the flame heat flux (q.../4, 


= 1) for the three thermoplastics tested. Similarly , Figures 6.2 and 6.3 show the plots 


for values of q.,,/G, equal to 2 and 3, respectively. It can be seen from these plots that 


there are only small variations in the solutions for the different materials for a given value 
of the ratio of external irradiance to flame heat flux. Figure 6.4 shows a plot of 
dimensionless mass loss rate versus dimensionless time after ignition for Nylon 6/6 at a 
number of ratios of external irradiance to flame heat flux. Figure 6.4 shows that the 
variation in the dimensionless mass loss rate for a given material is dependent on the 
values of the surface heat flux before and after ignition. However, it can be seen from 
Figures 6.1 - 6.4 that about 95% of the steady state mass loss rate is reached when the 
dimensionless time is approximately 0.3 - 0.5. The dimensionless time to ignition is 


represented by Equations (6.9) and (6.10), such that, 


c..,2 
1, = eed (6.12) 


The dimensionless mass loss rate at ignition, which is the y-intercept in Figures 


6.1- 6.4, can be represented as, 
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1 ] 
go | cee | aero aL (6.13) 
; r erica 


Examination of the dimensionless equations shows that, 


a AH, i, Tf, h(T,,-T,) 
Oe = (tty ot) 
q, Gs ds 
where, 
piste Lolly, h(T,, — To) 
L ,—, and -—~—*— 
ds dr 


are constant for a given material burning in the Cone Calorimeter. 

Figures 6.5 - 6.6 show plots of the calculated dimensionless mass loss rate as a 
function of dimensionless time along with the experimental data for Nylon. Similarly, 
plots for Polyethylene are shown in Figures 6.7 - 6.8, plots for Polypropylene are 
shown in Figures 6.9 - 6.10, and plots for PMMA are shown in Figures 6.11 - 6.12. 
These plots represent an alternative way of viewing the results. The scatter in the data in 


Figure 6.7 is attributed to bad data from the experiment. 
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Figure 6.1 Dimensionless mass loss rate as a function of dimensionless 
ume for q".,,/Q", = 1. 
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Figure 6.2 Dimensionless mass loss rate as a function of dimensionless 
time for q”.,,/Q"_ = 2. 
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Figure 6.3 Dimensionless mass loss rate as a function of dimensionless 
time for q"..,/q", = 3. 
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Figure 6.4 Dimensionless mass loss rate as a function of dimensionless 
time for Nylon. 
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Figure 6.5 Theoretical and experimental dimensionless mass loss rate as a function of 
dimensionless time for Nylon with a 35 kW/m? external irradiance. 
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Figure 6.6 Theoretical and experimental dimensionless mass loss rate as a function of 
dimensionless time for Nylon with an 80 kW/m? external irradiance. 
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Figure 6.7 Theoretical and experimental dimensionless mass loss rate as a function of 
dimensionless time for Polyethylene with a 36 kW/m? external irradiance. 
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Figure 6.8 Theoretical and experimental dimensionless mass loss rate as a function of 
dimensionless time for Polyethylene with an 87 kW/m? external irradiance. 
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Figure 6.9 Theoretical and experimental dimensionless mass loss rate as a function of 
dimensionless time for Polypropylene with a 27 kW/m? external irradiance. 
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Figure 6.10 Theoretical and experimental dimensionless mass loss rate as a function of 
dimensionless time for Polypropylene with a 65 kW/m? external irradiance. 
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Figure 6.11 Theoretical and experimental dimensionless mass loss rate as a function of 
dimensionless time for black PMMA with a 25 kW/m2 external irradiance. 


84 


PMMA 
75 kW/mA2 


m" L/q" 


S Experiment 
“== Calculated 


Figure 6.12 Theoretical and experimental dimensionless mass loss rate as a function of 
dimensionless time for black PMMA with a 75 kW/m? external irradiance. 
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CHAPTER 7 WOOD EXPERIMENTAL SETUP, 
PROCEDURE, AND MEASUREMENTS 


The wood experiments were performed in manner consistent with the testing of 
the thermoplastic materials described in chapter 3. The analysis of the wood results and 
the development of a model for charring materials 1s the focus of future work on this 
project. For now, only a description of the testing procedure and the results is 
presented. Further analysis of the results for both the mass loss rate and the heat flux is 
needed. Observations were taken during the experiments and were recorded in a 


notebook to aid in the analysis. 


7.1 Experimental System 

Ignition and burning rate experiments were performed using a radiant cone heater 
assembly as shown in Figure 3.1. Thermocouple, heat flux meter, and mass loss 
readings were recorded using a LabView data acquisition system and a Macintosh IIFX 
computer. Data was recorded every second throughout the experiments. The data 
acquisition system was compiled using the simple virtual instruments (VI’s) built into 


LabView, however, a copy of the program is not included in this report. 


7.2 Wood Samples 

Two materials were selected for analysis in this study, namely; Red Oak and 
Redwood. The Redwood samples measured 10 cm (4 in.) x 10 cm (4 in.) x 1.91 cm 
(0.75 in.) thick. The Red Oak samples measured 10 cm (4 in.) x 10 cm (4 in.) x 3.175 
cm (1.25 in.) thick. The samples were conditioned in an atmosphere of approximately 
20 °C and 50% relative humidity for at least two weeks prior to testing. Aluminum foil 
was wrapped around the sides of the samples in an attempt to inhibit edge burning. This 


did not completely eliminate edge burning effects, but it did significantly reduce the 
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effects of edge burning to maintain one dimensional burning behavior. 
Samples were placed on a standard metal cone holder in the horizontal orientation 
on a bed of Kaowool so that heat loss effects through the back of the samples remained 


as insignificant as possible. 


7.3 Experimental Procedure 

The experimental procedure consisted of exposing a sample, in the horizontal 
orientation, to a constant external irradiance from the cone heater assembly. The time to 
piloted ignition was measured and mass loss data recorded for each test. Prior to the 
onset of each experiment the cone heater needed to be set to the desired constant external 
irradiance. The initial incident heat flux was determined using a Medtherm 2.5 cm (1 
inch) diameter thermopile type heat flux transducer (Model # 64-5SB-20) situated such 
that it was in the same location that the center of the sample would be, and so that the 
face of the gauge was 2.5 cm (1 inch) below the base of the cone heater. The experiment 
was not started until a constant heat flux recording was obtained for at least a one minute 
period. The et flux gauge could not be left in place during the experiment. 
Subsequently, the constancy of the cone heater was ensured by monitoring the 
temperature of the cone coil. Two thermocouples were located in the cone coil. The 
thermocouple readings were averaged and displayed every second so that the cone heat 
flux, which corresponds to a given coil temperature, could be monitored. The cone coil 
temperature was kept constant (+ 5° C) by manually adjusting the current to the coil 
using a 220 Volt transformer. This approach allowed for adjustment of the current to 


compensate for temperature fluctuations induced by flame impingement on the cone coil. 
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7.3.1 Ignition 


The arrangement of the assembly was such that the top of the sample was initially 
located 2.5 cm (1 inch) below the base of the cone heater as is done for standard Cone 
Calorimeter tests [3]. The load cell was oriented between two guide bars to ensure 
proper placement of the sample underneath the center of the cone heater assembly. A one 
inch methane flame, located on one edge approximately 1.25 cm (1/2 inch) above the 
surface of the sample, was used as a pilot ignition source. Figure 3.3 shows the 
arrangement of the pilot and the assembly. 

Once the cone was set to the desired constant external irradiance and the sample 
was prepared and situated on the load cell, the load cell was shifted into position 
underneath the cone. Simultaneously, the data acquisition system was initiated. It 
should be noted that sliding the sample into place caused some bouncing of the load cell. 
This effects the mass loss recordings for the first few seconds of the experiment. This is 
not however expected to have adverse effects on the results since ignition nominally took 
longer than ten seconds. 

The ignition time is defined as the time at which a continuous flame is supported 
on the material surface. In some instances flashing occurred on the surface of the sample 
prior sustained flaming. However, in all cases the ignition time was taken as the time at 
which flaming was sustained over the entire surface of the specimen. Ignition times for 
all of the wood tests for irradiances ranging from 20 kW/m2 to about 80 kW/m2 are 
shown in Appendix G. 

In some experiments a fine wire Type K thermocouple (0.003 in. diameter) was 
mounted on the thermoplastic material surface to attempt to measure the surface 
temperature at ignition. The-placement of the thermocouples on the wood samples was 
done in a manner consistent with a method recommended by Atreya [18]. The wires and 


the bead were flattened to obtain a film thermocouple approximately 0.001 inches thick. 
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Figure 7.1 Surface temperature of Redwood as a function of time with a 30 kW/m2 
external irradiance. 
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A very fine incision was made in the surface of the sample and the thermocouple was slid 
underneath the skin of the incision which was approximately 0.001 inches thick. A 
small amount of wood glue was used to secure the rest of the thermocouple. The entire 
assembly was then pressed and allowed to set for approximately one hour. 
Measurements of the surface temperature were recorded until the onset of burning. 
Figure 7.1 shows a typical result of the measured surface temperature for the preheating 
period for Redwood exposed to a 30 kW/m? external irradiance. The temperature at the 
onset of the sudden rise in temperature is defined as the ignition temperature (Tj,). 
Results for the surface temperature measurements for the wood experiments at various 


external irradiances are shown in Appendix H. 


13.2 Mass Loss Rate 
Experiments to determine the mass loss rate of the materials were performed 
concurrently with the ignition experiments. The mass loss of the samples was recorded 
using a load cell. The mass loss readings were recorded every second for approximately 
a 1200 second period. 
The transient mass loss rate is found using a seven point least squares fit of the 
mass loss data. The mass loss rate at some time t; is given as, 
i+ 2 i+ 2 i+ 2 
7 ny (mt, 7 y (m,) » (t) (7.1) 
n= i-2 n=i-2 n= i-2 
1'4-22 
1 hy Se ) = 


n=i-2 


where, 1 is a given measured data point, 
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Figure 7.2 Transient mass loss rate of Redwood with a 36 kW/m? external irradiance. 
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m is the mass at 1, 
and, t is the time at i. 

Figure 7.2 shows a typical mass loss rate result for Redwood exposed to an 36 
kW/m? irradiance. Results for Redwood at other exposures and for Red Oak are shown 
in Appendixes I and J respectively. It should be noted that a seven point least squares fit 
was used for the calculation of the mass loss rates for the wood experiments whereas a 
five point least squares fit was used for the thermoplastic experiments. The reason for 
the increased fit was due to scatter in the data. A low pass filter consisting of a 10 kQ 
resistor and an 8 uF capacitor was used to help reduce some of the noise. Data points 
were taken every second for the experiments, however, the recommended approach for 
future experiments using LabView is to take a number of data points, say 100, per 
second and to have LabView average these points and return a single data point at each 
second. This would yield much better results and is recommended for future 
experiments. In addition, use of the intermediate and advanced VI’s built into LabView 
would allow for filtering of the data to be done with the LabView system. 

In addition to the use of the seven point least squares fit, it was necessary to 
perform a five point smooth on the mass loss data. Because of this, the peak mass loss 
rates for the wood experiments were determined manually by plotting the mass loss 
versus time for each experiment and taking the slope at the time the maximum occurred. 
This enabled for a more accurate determination of the peak mass loss rate. The plots in 


Appendices I and J contain the corrected peak mass loss rates. 


7.3.3 Flame Heat Flux 


The flame plus external heat flux was measured during some of the experiments 
using a 3.175 mm (1/8 inch) Medtherm thermopile type heat flux gage (Model # 8-1-10- 


4-0-36-20680K). Due to contraction of the wood against the gage causing erroneous 
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Figure 7.3 Total incident heat flux, measured by sensor, of Red Oak with a 28 
kW/m2 external irradiance. 


93 


mass loss measurements, it was not possible to measure the mass loss and the heat flux 
to the surface in the same experiments. Subsequently, separate tests at approximately the 
same external irradiance were performed to measure the mass loss and the heat flux 
respectively. To avoid excessive condensation and re-evaporation of pyrolisates on the 
surface of the heat flux gauge the water temperature was raised to approximately 65 °C. 
This was accomplished by circulating water through a copper coil which was heated by a 
methane flame. While raising the temperature of the water to 65 °C may have reduced 
the amount of condensation on the sensor, some slight condensation and re-evaporation 
did occur. The gauge was insulated on the sides with Kaowool to help minimize 
extraneous results. Figure 7.3 shows a typical result for the heat flux measurement for 
Red Oak exposed to a 28 kW/mz2 external irradiance. Results for the other experiments 
where the heat flux was measured are shown in Appendix K for Redwood and Appendix 
L for Red Oak. 


7.4 Calibration 
Calibration of the load cell, thermocouples, and heat flux gages was necessary 
before any experiments were conducted and only occasionally once testing began. The 


methods used to calibrate each of the instruments are discussed below. 


1.4.1 Load Cell 
The load cell was calibrated by first adjusting the weight of the sample holder and 


insulation to approximately zero. At this point, standard weights between 1 and 100 
gram were added to the load cell and the voltage readings were recorded. The 
relationship between weight and voltage was found to be linear, with a calibration 
constant of 582.0 g/V. Since only mass loss rate is desired, only the difference in mass 


is necessary. Therefore, small fluctuations in the zero weight did not affect the results. 
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At the start of each day, the load cell was rechecked for accuracy using 1 and 2 gram 
weights. Furthermore, in some instances the mass of the sample was measured using a 
scale and the results compared to that indicated by the load cell. No significant errors 


were found for the mass loss measurements. 


7.4.2 Heat Flux Gages 
The 2.5 cm (1 inch) and 3.175 mm (1/8 inch) heat flux gages were initially 


calibrated using the standard NIST heat lamp calibration box. They were tested with 
water at room temperature and again with water at 46 °C. No notable differences in the 
calibration constant were observed. Although the water temperature did not affect the 
calibration constants, it did affect the zero reading. Unlike the mass loss, the heat flux 
measurement must be zero when there is no incident heat flux. The hot water resulted in 
an initial negative heat flux measurement due to convection between the sensor and the 
atmosphere. The heat flux measurement was calibrated to be zero, under no external 
irradiance, before each experiment in order to obtain consistent results. Besides the zero 
flux Pea eeTent the 2.5 cm (1 inch) gage did not require calibration before each of the 
éxperiments because it was positioned such that it did not get disturbed or touched. 
However, the 3.175 mm (1/8 inch) heat flux meter often became dirty after an 
experiment. Subsequently, it was calibrated before each experiment using the 2.5 cm (1 
inch) gage as a standard. Differences between the 3.175 mm (1/8 inch) heat flux gage’s 
measurement before and after each experiment were found to be within 5 percent. The 
calibration constant for the 2.5 cm (1 inch) heat flux meter was determined to be 5235 


kW/m2 V and the calibration constant for the 3.175 mm (1/8 inch) heat flux gage was 


Aeremminedito be 1184. kWwinaive 
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7.4.3 Thermocouples 


The thermocouples used to measure the coil temperature and surface temperature 
were calibrated before any testing began. This was done by nace them in an ice bath 
and also in boiling water for a prolonged period of time. Temperature readings were 
found to be accurate within + 2 °C. The thermocouples in the coil did not require 
additional calibration as the absolute temperature was not necessary. The thermocouples 
for the surface were calibrated against each other before each experiment. 

Although the experiments yielded extraneous results for some tests, no major 
calibration errors were found. Inconsistencies in the results are likely to be due to 


factors other than calibration errors. 
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CHAPTER 8 CONCLUSIONS 


The modeling prescription developed by Quintiere and Rhodes [8], which can 
utilize Cone Calorimeter data to derive useful properties needed to predict ignition and 
transient burning rates for thermoplastic-like materials, has been shown to yield good 
results. The level of accuracy has been demonstrated for Nylon 6/6, Polyethylene, 
Polypropylene, and PMMA. The simplicity of the ignition model is advantageous and 
can be used to infer the critical flux for ignition without direct measurement. 

Heat flux to the material surface was not measured for the present experiments, 
rather the flame heat flux to the surface was inferred using the method described in 
section 5.2. The inference of a constant flame heat flux to the surface for thermoplastic- 
like materials burning in the Cone Calorimeter has been shown to yield good results. 
Values for the flame heat flux were found to be 30 kW/m? for Nylon, 25 kW/m2 for 
Polyethylene, 14 kW/m? for Polypropylene, and 37 kW/m? for PMMA. The constancy 
of the flame heat flux appears to be an attribute of all thermoplastics burning in the Cone 
Calorimeter due to the long column like shape of the flame. 

For all of the thermoplastic materials examined thus far, the calculated ignition 
temperature has been found to be less than the measured ignition temperature. This is 
expected to be an attribute of the modeling prescription and the results are not expected to 
be void of physical significance. Use of the calculated ignition temperatures in the model 
has been shown to yield good results. | 

Ignition, burning rate, and heat flux measurements have been taken for Redwood 
and Red Oak burning in the Cone Calorimeter. An assessment of the results and an 
attempt to formulate a modeling prescription for charring materials, similar to the 
prescription for thermoplastics, should be attempted. The intent is to maintain simplicity 


while maintaining completeness of the important processes. The framework for the 
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charring model can be found in the formulation by Quintiere [1]. 
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APPENDIX A THERMOPLASTIC DATA ACQUISITION PROGRAM 
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10 OPTION BASE 1 

20 DIM In$[100}, Out(6) 

Zz =-5.115307103E-2 

22 P1=2.485028007E+4 

23. = P2=-3.826622822E+5 

24 P3=9.966105673E+7 

25 P4=-1.082062357E+10 

26 P5=6.039285524E+11 

uF | P6=-1.910899962E+13 

28 P7=3.478234730E+14 

29 P8=-3.399102821E+15 

30 P9=1.382851398E+16 

cy, CALL Set_datacom 

33 CLEAR 722 

34 CLEAR 709 

35 FOR I=1 TO 8 

36 ON KEY I LABEL “ ““GOSUB Beeper 

37 NEXT I 

38 OUTPUT 709;’F00 L05 11” 

39 ~ OUTPUT 709;700” 

40 OUTPUT 722;”L1 RS1 DO ZO F1 R1 0.01STD SO1 10.0STI 6STN T4 SM004 - 
7STR RER Q” 

41 OUTPUT 722;’X1” 

42 Oncom=0 

43.  ONKEY 1LABEL “COM OFF“ GOSUB Comonoff 

44 ON KEY 4 LABEL “END TEST “ GOTO Endit 

45 ON KEY 5 LABEL “START TEST “ GOTO Onward 

46 ! 

47 ! 

48 Spin: ! 

49 GOTO Spin 

50 ! 

51 ! 

52 Onward: Ti0=TIMEDATE MOD 86400 

53 Restart: ! 

54 TRIGGER 722 

55 TI1=(TIMEDATE MOD 86400)-TI0 

58 ENTER 722;In$ 

59 OUTPUT 722;”X1” 

60 ENTER In$;Out(*) 

61 IF Oncom=1 THEN 

62 FOR I=2 TO 6 STEP 2 

63 Cd)=Out(1)*(P6+Out(I)*(P7+Outr()*(P8+Out(1)*P9))) 

64 TdM)=P0+Out(l)*(P1+Out(I)*(P2+Out()*(P3+Out(1)*(P4+Out(I)*(P5+C())))) 

65 NEXT I 

66 OUTPUT 9 USING “DDDD.D,?”;Til 

69 OUTPUT 9;” “; 

70 OUTPUT 9 USING ‘“DD.D,#’;Out(1)*5465 

ll OUTPUT 9;” 

104 ON ERROR GOTO Error 
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OUTPUT 9 USING “DDDD._D,#”;T(2)-2.3 


100 


ip GUIPULS:7 . *: 

76 OUTPUT 9 USING “DDD.D,#”;Out(3)* 11360 

ry OUirur 3 

78 OUTPUT 9 USING “DDD.D,#”;T(4) 

79 OUTPUT 9;” “; 

80 OUTPUT 9 USING “DDDD.DD,#”’;Out(5)* 1000*(-0.58597) 
81 OUTPUT 9;” “; 

82 OUTPUT 9 USING “DDDD.D”;T(6) +0.4 

83 ENDIF 

84 PRINT USING “DDDD.D,#’;Til 

85 PRINT “ “;Out(*) 

86 LOOP 

87 EXIT IF (TIMEDATE MOD 86400)-Ti0>=Til +2 

88 END LOOP 

89 GOTO Restart 

90 Comonoff: ! 

91 IF Oncom=0 THEN 

92 Oncom=1 

93 ON KEY 1 LABEL “COM ON “ GOSUB Comonoff 
94 ELSE 

95 Oncom=0 

96 ON KEY 1 LABEL “COM OFF “ GOSUB Comonoff 


97 ENDIF 

98 RETURN 

99 Beeper: ! 

100 BEEP 

101 RETURN 

102 Error: ! 

103 a. QUTRUT DO: sites 
104 GOTO75 

106 Endit: ! 

107 ABORT7 


110 END 

Ap Lie! 

112 SUB Set_datacom 
113, ~Sc=9 


114 CONTROL Sc,0;1 
115 CONTROL Sc,3;9600 
124 CONTROL Sc,4;3 
134. CONTROL Sc,5;0 
144 SUBEND 
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APPENDIX B IGNITION TIME FOR THERMOPLASTICS 
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Test Number External Irradiance Time to Ignite 
(kW/m2) 


Table B1 Experimental ignition time for Nylon 6/6 
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Table B2 Experimental ignition time for polyethylene. 
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Table B3 Experimental ignition time for polypropylene. 
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APPENDIX C MEASURED SURFACE TEMPERATURE 
FOR THERMOPLASTICS 


106 


| Mater Material/Test Number | External | External Irradiance | Measured Ignition Te Temperature 
| (kW/mz2) 
Nylon/n8 31 500 


[Poiyetylenees | 9ST SSC 

oe 

oS {aes 
Polypropylene/pp9 34 315 


TABLE C1 Measured ignition temperature for thermoplastics. 
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Polyethylene - 59 kW/mA2 
Ignition at 35 sec. 
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Figure C1 Surface temperature results for Polyethylene with a 59 kW/mz2 external 
irradiance. 
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Polypropylene - 61 kW/m42 
Ignition at 22 sec 
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Figure C2 Surface temperature results for Polypropylene with a 61 kW/m2 external 
irradiance. 
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Polypropylene - 34 kW/m42 
Ignition at 66 sec. 
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Figure C3 Surface temperature results for Polypropylene with a 34 kW/m2 external 
irradiance. 
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Nylon - 31 kW/m42 
Ignition at 505 sec 
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Figure C4 Surface temperature results for Nylon with a 31 kW/m? external 
irradiance. 
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APPENDIX D NYLON MASS LOSS RATE 
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[ Test Number Number | External Irradiance | lance | Steady "Steady State Mass Loss Rate | Mass Loss Rate 
— oe Ss) 
n2 


x 


n3 Pree 25 


TABLE D1 Steady state mass loss rate for Nylon 6/6. 
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Figure D1 _ Transient mass loss rate results for Nylon with a 31 kW/mz2 external 
irradiance. 
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Nylon - 31 kW/mA2 
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Figure D2 Transient mass loss rate results for Nylon with a 31 kW/m2 external 
irradiance. 
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Nylon - 35 kW/m42 
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Figure D3 __ Transient mass loss rate results for Nylon with a 35 kW/m? external 
irradiance. 
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Nylon - 50 kW/m42 
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Figure D4___ Transient mass loss rate results for Nylon with a 50 kW/m external 
irradiance. 
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Nylon - 60 kW/m42 
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Figure D5 __ Transient mass loss rate results for Nylon with a 60 kW/m? external 
irradiance. 
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Nylon - 75 kW/m42 
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Figure D6 Transient mass loss rate results for Nylon with a 75 kW/m? external 
irradiance. 
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APPENDIX E POLYETHYLENE MASS LOSS RATE 
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| Test Number Steady State Mass Loss 
: (g/m2 s) 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 


TABLE E1 Steady state mass loss rate for Polyethylene. 
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Polyethylene - 26 kW/m42 
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Figure El Transient mass loss rate results for Polyethylene with a 26 kW/m? external 
irradiance. 
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Polyethylene - 46 kW/mA2 
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Figure E2 Transient mass loss rate results for Polyethylene with a 46 kW/m? extemal 


imradiance. 
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Polyethylene - 59 kW/m42 


Mass Loss Rate (g/m‘“2-s) 


0 200 400 600 800 1000 1200 


Time (sec) 


Figure E3 Transient mass loss rate results for Polyethylene with a 59 kW/mz2 external 
irradiance. 
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Polyethylene - 61 kW/m42 
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Figure E4 Transient mass loss rate results for Polyethylene with a 61 kW/m? external 
irradiance. 
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Polyethylene - 70 kW/m42 
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Figure E5 Transient mass loss rate results for Polyethylene with a 70 kW/m2 external 
irradiance. 
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Polyethylene - 87 kW/m42 
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Figure E6 Transient mass loss rate results for Polyethylene with a 87 kW/m? external 
irradiance. 
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APPENDIX F POLYPROPYLENE MASS LOSS RATE 
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Test } Number External | External Inradiance | Steady State Mass Loss Rate State Mass Loss Rate 
“se eee eer S) 


TABLE F1 Steady state mass loss rate for Polypropylene. 
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Mass Loss Rate (g/m‘“2-s) 


Figure Fl 


Polypropylene - 20 kW/mA2 


— Cpt. = 

acligze re 
oO. 
t- 


0 200 400 600 800 1000 1200 


Time (sec) 


Transient mass loss rate results for Polypropylene with a 20 kW/m2 external 
irradiance. 
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Transient mass loss rate results for Polypropylene with a 27 kW/m? external 
irradiance. 
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Polypropylene - 34 kW/m42 
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Figure F3 Transient mass loss rate results for Polypropylene with a 34 kW/m2 external 
irradiance. 
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Polypropylene - 39 kW/m42 
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Figure F4 Transient mass loss rate results for Polypropylene with a 39 kW/m? external 
irradiance. 
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Polypropylene - 50 kW/m42 
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Figure F5 Transient mass loss rate results for Polypropylene with a 50 kW/mz2 external 
irradiance. 
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Polypropylene - 61 kW/m42 
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Figure F6 Transient mass loss rate results for Polypropylene with a 61 kW/m2 external 
irradiance. 
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Polypropylene - 61 kW/m42 
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Figure F7 Transient mass loss rate results for Polypropylene with a 61 kW/m? external 
irradiance. 
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Polypropylene - 65 kW/m42 
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Figure F8 Transient mass loss rate results for Polypropylene with a 65 kW/m? external 
irradiance. 
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Polypropylene - 70 kW/m42 
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Figure F9 Transient mass loss rate results for Polypropylene with a 70 kW/mz2 external 
irradiance. 
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APPENDIX G IGNITION TIME FOR WOODS 
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g|& 
=| < 
duh 


{ 
= 
Nw 


Maverial Test Number 
| 

| 

| 

| 

| 

| 

| 

| Redwood/Rw7 
ee 
| 

| 

| 

| 3 


Red ne 


TABLE G1 Experimental ignition times for Redwood and Redoak. 
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APPENDIX H MEASURED SURFACE TEMPERATURE FOR WOODS 
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[  Material/Tes est t Number | External Irradiance | External Irradiance || Measured Ign Measured Ignition 
| KW?) dee ie _ Temperature G Oem 


eee 3) | 


Redwood/RW 14 | 


Redwood/RW 15 


Redwood/RW 16 


Table H1 Measured ignition temperature for Redwood and Red Oak. 
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Figure H1 Surface temperature results for Redwood with a 21 kW/m? external 
irradiance. 
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Figure H2 Surface temperature results for Redwood with a 42 kW/m? external 
irradiance. 
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Figure H3 Surface temperature results for Redwood with a 54 kW/m? external 
irradiance. 
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Figure H4 Surface temperature results for Red Oak with a 27 kW/m2 external 
irradiance. 
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Red Oak 
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Figure H5 Surface temperature results for Red Oak with a 52 kW/m? external 
irradiance. 
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Red Oak 
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Figure H6 Surface temperature results for Red Oak with a 76 kW/m? external 
irradiance. 
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APPENDIX I REDWOOD MASS LOSS RATE 
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Figure I1 = Mass loss rate results of Redwood with a 20 kW/m? external irradiance. 
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Redwood 
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Figure I2 Mass loss rate results of Redwood with a 21 kW/m? external irradiance. 
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Figure I3_ Mass loss rate results of Redwood with a 30 kW/m2 external irradiance. 
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Figure I4 Mass loss rate results of Redwood with a 36 kW/m external irradiance. 
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Mass loss rate results of Redwood with a 42 kW/m2 external irradiance. 


Figure I5 
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Figure I6 Mass loss rate results of Redwood with a 54 kW/m2 external irradiance. 
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Redwood 
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Figure I7 _— Mass loss rate results of Redwood with a 55 kW/m? external irradiance. 
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Figure I8 = Mass loss rate results of Redwood with a 60 kW/m? external irradiance. 
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APPENDIX J RED OAK MASS LOSS RATE 
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Figure Jl Mass loss rate results of Red Oak with a 27 kW/m? external irradiance. 
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Mass loss rate results of Red Oak with a. 46 kW/m2 external irradiance. 


Figure J2 


160 


Red Oak 
76 kW/m42 


Mass Loss Rate (g/m42-s) 


0 200 400 600 800 1000 1200 


Time (sec) 


Figure J3 Mass loss rate results of Red Oak with a 76 kW/m? external irradiance. 
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APPENDIX K REDWOOD INCIDENT HEAT FLUX 
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Figure K1 Incident heat flux results for Redwood with a 20 kW/m? external irradiance. 
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Figure K2 Incident heat flux results for Redwood with a 23 kW/m2 external irradiance. 
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Figure K3__ Incident heat flux results for Redwood with a 30 kW/m? extemal irradiance. 
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Figure K4 __ Incident heat flux results for Redwood with a 42 kW/m? external irradiance. 
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Figure K5__ Incident heat flux results for Redwood with a 44 kW/m? external irradiance. 
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Figure K6 ‘Incident heat flux results for Redwood with a 55 kW/m2 external irradiance. 
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Figure K7 Incident heat flux results for Redwood with a 76 kW/m? external irradiance. 
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APPENDIX L RED OAK INCIDENT HEAT FLUX 


170 


Red Oak 
28 kW/mA2 


Flux (kW/m*2) 


0 200 400 600 800 1000 1200 


Time (sec) 


Figure L1 _ Incident heat flux results for Red Oak with a 28 kW/m external irradiance. 
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Figure L2 Incident heat flux results for Red Oak with a 49 kW/m? external irradiance. 
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Figure L3 Incident heat flux results for Red Oak with a 80 kW/m? external irradiance. 
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NOMENCLATURE 


c - specific heat 

h, - convective heat transfer coefficient 
k - thermal conductivity 

1 - beam length 

L - heat of gasification 

m - mass 

Q - power output 

q - heat flow 

I - stoichiometric oxygen to fuel mass ratio 
T - temperature 

t - time 

y - Space coordinate 

y-dimensionless mass loss rate (ch. 6) 


Yox, « - ambient oxygen mass fraction 
@ - thermal diffusivity 

6 - thermal penetration depth 

AH, - heat of vaporization 

AH, - heat of combustion 

€ - emissivity 

k - absorption-emission coefficient 


p - density 
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6 - Stefan Boltzmann constant 
X - heat fraction 

€ — variable (Eq. 4.16) 

t — dimensionless time 


 - variable (Eq. 6.2) 


x - variable (Eq. 6.3) 
Subscripts 

Cc - convective 

cr - critical 

ext - external 

fl - flame 

fl,c - flame convection 
fl,r - flame radiation 
§ - gas 

ig - ignition 

O - initial, ambient 

m - mean 

r - radiative 

S - steady 

V - vaporization 
Superscripts 

( ) - per unit time 


( )” - per unit area 
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